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Abstract 

We consider the valence-bond-solid ground state of the q-deformed higher-spin 
AKLT model (g-VBS state) with q real. We investigate the eigenvalues and eigen- 
vectors of a matrix (G matrix), which is constructed from the matrix product rep- 
resentation of the g-VBS state. We compute the longitudinal and transverse spin- 
spin correlation functions, and determine the correlation amplitudes and correlation 
lengths. 

1 Introduction 

In one-dimensional quantum systems, a completely different behavior for the integer spin 
chains from the half-integer spin chains was predicted by Haldane [TJ [2] . The antiferro- 
magnetic isotropic spin-1 model introduced by Affleck, Kennedy, Lieb and Tasaki (AKLT 
model) [3], whose ground state can be exactly calculated, has been a useful toy model to 
validate Haldane's prediction of the massive behavior for integer spin chains. Moreover, it 
lead to a deeper understanding for integer spin chains such as the discovery of the special 
type of long-range order [HE]. The AKLT model has been generalized to higher-spin mod- 
els, anisotropic models, etc |Sl El El El EH HH H21 H31 H31 HSl HEl - Th e Hamiltonians are 
essentially linear combinations of projection operators with nonnegative coefficients, and 
their ground states are called valence-bond-solid (VBS) state. Recently, the VBS state 
was investigated in perspective of its relation to quantum information and experimental 
implementation by means of optical lattices, see Refs. |18[ [19] for example. 

There are largely three types of representations for the ground state which are equiv- 
alent to each other: the Schwinger boson representation, the spin coherent representation 
and the matrix product representation. For isotropic higher-spin models, the spin-spin 
correlation functions [20] and the entanglement entropy [21[ [22] have been calculated by 
utilizing the spin coherent representation and the properties of Legendre polynomials. For 
the g-deformed spin-1 model, spin-spin correlation functions were evaluated [7] [8] [9] from 
the matrix product representation. 

In this paper, we consider the ground state of a g-deformed higher-integer-spin model 
which was constructed recently in Ref. |24] (g-VBS state). From its matrix product repre- 
sentation, we analyze one and two point functions of the g-VBS ground state. We notice 
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that a matrix, which is constructed from the matrix product representation, plays a funda- 
mental role in computing correlation functions, especially spin-spin two point correlation 
functions. Investigating the structure of the matrix in detail, we obtain its eigenvalues and 
eigenvectors. Utilizing the results, we determine the correlation amplitudes and correlation 
lengths of the longitudinal and transverse spin-spin correlation functions. 

This paper is organized as follows. In the next section, we briefly review the quantum 
group U q (su(2)), and investigate the finite dimensional highest weight representation in 
terms of Schwinger bosons. In Section [31 we precisely define the higher-spin generalization 
of the g-deformed AKLT model on an L-site chain, and rigorously derive its g-VBS ground 
state in a matrix product form. The squared norm of the state will be written in terms of 
the trace of the L-th power of a matrix G, which plays an important role in this paper. In 
section HI we obtain the eigenvalues and eigenvectors of G. Utilizing them, we compute one 
and two point functions in Section [5j Especially, we determine the correlation amplitudes 
and correlation lengths of the longitudinal and transverse spin-spin correlation functions. 
Section [6] is devoted to the conclusion of this paper. 



2 The quantum group U q (su(2)) 



We introduce several notations, fixing a real number q throughout this paper. Let us 
define the q- integer, q- factorial and g-binomial coefficient for TV" € Z>o as 



[N] 



q N -q- N 
q-q 1 



N 
K 



[N]\ 
[N]\ 



f N 

i=i 

1 N = 0, 

K = 0,...,N, 



(2.1) 



[K]\[N -K]\ 

otherwise, 



respectively. 

The quantum group U q (su(2)) [25j ES] is defined by generators X + ,X~ and H with 
relations 



[x + ,x- 



q-q- 



[H.X^] = ±2X d 



(2.2) 



The comultiplication is given by 

A(X ± ) = X ± ®q H/2 + q- H/2 (S)X ± , A(H) = H ® Id + Id ® H. (2.3) 
U q {su{2)) has the Schwinger boson representation, where the generators are realized 



as 



X + = a ] b, X- = b ] a, H = N a — N b , 

with g-bosons a and b satisfying 

aa j - qa j a = q~ Na , bb j - qb j b = q~ Nb , 

[N a ,a] = -a, [N a ,at} = a\ [N b ,b} = -b, [JV 6s 6+] = &+. 



(2.4) 

(2.5) 
(2.6) 
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We denote the space where (2j + l)-dimensional highest weight representation of U q (su(2)) 
is realized by Vj. The basis of Vj is given by 



\j;m) 



(at)i+ m (fet)i-m 
+ m]\\j -m}\ 



\vac), {m = 



(2.7) 



The Weyl representation which we describe below, is an equivalent representation to 
the Schwinger boson representation, and is efficient for practical calculation. Let us denote 
the g-bosons a and b acting on the ct-th site as a a and b a . The Weyl representation is to 
represent a a , ba, a a and b a on the space of polynomials C[x a ,y a ] as 



O'a — ^cn &a — Von 



1 U Q U n -1 



x a q-q 



.1 h 

_1 J u a 



1 D y q " - Dp, 

y a q- q~ l 



where 



Dp a f(x a ,y a ) = f(px a ,y a ), Df* f(x a ,y a ) = f(x a ,py a ). 
The generators of U q {su{2)) are now represented as 



X+ 



y a g-q 



_i > x a 



Vat ^1 



X (-. 



q 



D x n a D ya x . 

1 q 1 



(2. 



(2.9) 



(2.10) 



The tensor product of two irreducible representations has the Clebsch-Gordan decom- 
position 



25' 



Vs®V s = Q)Vj 
j=o 

2S 

\S;mi) <g) \S;m 2 ) = } j 



S S J 

mi m 2 mi + m 2 



\J;mi + m 2 ), 



(2.11) 
(2.12) 



where 



Si S 2 J 
mi m 2 m 



/_i\Si-mi m 1 (m 1 +m 2 +l)+{S 2 (S2+l)-S , i(5i+l)-J(J+l)}/2 
,m\ L ) q 



[J + m]\[J - m]\[Si - mi]\[S 2 - m 2 ]![5i + S 2 - J]! [2 J + 1] 
[Si + mi]\[S 2 + m 2 ]![5i - S 2 + J]\[S 2 - Si + J]! [Si + S 2 + J + 1]! 



(2.13) 



Min( J— m,Si— mi, 52+ J— mi) 



2=Max(0, — Si— mi, J— S2— mi) 



(-q m+J+1 ) z [Si + mi + z]l[S 2 + J - mi - z}\ 
[z}\[J - m - z]\[Si - mi - z]l{S 2 - J + mi + z]V 



is the ^-analog of the Clebsch-Gordan coefficient |23| . (The factor g m i m 2/ 2 j s missing in 
Ref. [23J.) This coefficient is compatible with the inverse of the decomposition (I2.12p 



\J;m)= 



mi+rri2=m 



S S J 
mi m 2 mi + m 2 



15; mi) ® \S;m 2 ). 



(2.14) 
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For later purpose, we will also investigate the Clebsch-Gordan decomposition of U q (su(2)) 
in terms of the Schwinger boson or the Weyl representation. Utilizing 



(2.15) 



one can show that the highest weight vector vj E Vj (AX + vj = 0) acting on the a-th 
and /3-th sites is given by 



2S-J 



v J = (x a x l3 ) J Y\ {xayp - q 2{y s l) Xf3y a ). 



(2.16) 



Moreover, we can show the following: 
Proposition 2.1. 



fi=0 



2fiS 


' J ' 




J 








n — fi 



2S-J 



(2.17) 



A proof of this proposition is given in Appendix IA"1 
Remark 2.2. Let n > 2J + 1. Noting 



" J " 




J 


. V- _ 




n — fj, 



0, 



for < n < n, one can see that 



(2.18) 



(2.19) 



3 g-VBS state 

The model we treat in this paper is an anisotropic integer spin-5 Hamiltonian on an L-site 
chain with the periodic boundary condition 



25' 



n = E E Cj(k,k + l)(7rj] 
k& L J=S+1 



k,k+l 



(3.1) 



where Cj(k, k + 1) > 0, and (vrj) fc k+1 , which acts on the k-th and (k + l)-th sites, is the 
Uq(su(2)) projection operator from Vs <S> Vs to Vj as 



S S J 

mi m 2 m i + ^2 



S S J 
m'i m' 2 m'i + m' 2 



vrj= 

mi ,rri2,mj ,m' 2 =— S 

x ^m 1 +m 2 ,m' 1 +m' 2 \S;m' l ){S]m 1 \ ® \S; m' 2 ) (S; m 2 |. 



(3.2) 
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The nonnegativity 

W ( 7r >/)fc,fc+i IV*) > ( for an y vector \ip)) (3.3) 

implies that all the eigenvalues of T~L are nonnegative. (Of course, is the Hermitian 
conjugate of |^).) Moreover, we will see that the energy of the ground state |^) is zero: 

H\V) = 0. (3.4) 

Since we set Cj(k, k + 1) > 0, we find that (|3.4p is equivalent to 

(vrj) fcife+1 |*) = (V& £ Z L , VJ € {S + 1, . . . , 25}) , (3.5) 

noting the nonnegativity (|3.3|) . From Proposition 12.14 one observes that any vector in 
©o< J<s Vj C V^f ® Vg of the fc-th and (fc + l)-th sites has the form 

s 

£ C AB xty S k - A ^k + iy S k +i H(q m x k y k+1 -q- m y k x k+1 ), (3.6) 

0<A,B<S m=l 

where Cab does not depend on x k ,y k , x k+ \ or y k +\. Thus, the condition (|3.5p imposes 
the restriction that \^f) has the form 

s 

|tf> = P({x k } ke z L ,{y k } keZL ) Yl Y[(q m x k y k+1 -q- m y k x k+1 ) (3.7) 

k£Z L m=l 

with some polynomials P such that this form is consistent with (|3.6|) for Vfc G Z^. The 
unique choice of P with such consistency is a constant (which can be set to be 1), and we 
achieve the unique ground state 

S 

l*>= II W^^kVk+i-q^ykXk+i). (3.8) 

fceZj, m=l 

In the Schwinger boson representation, we have 

l*)=IIII (« m 4&l +1 - q- m b\a\ +l )\^ C ), (3.9) 
fceZj, m=l 

which is a generalization of the q = 1 case [5J. Note that each site have the correct spin 
value: N k \V) = S\^) (k € Z L ) where N k := (N ak + N bk )/2. Our ground state is a g- 
deformation of the valence-bond-solid (VBS) state, which we call g-VBS state, see figure 

The Schwinger boson representation of the ground state (|3.9p can be transformed 
into the following equivalent form called the matrix product representation [23], which 
generalizes the q = 1 [27] or 5 = 1 [7J case. Noting (|2.7p . we have 

|*) =Tr[ 5l * 5 2*---*0L-i*PL], (3-10) 
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Figure 1: Conceptual figure of the g-VBS state. Each line is a g-deformed valence bond, 
and the circle Q represents the g-symmetrization of spin- 1/2 particles • at each site. 



where g k is an (5 + 1) x (5 + 1) vector- valued matrix acting on the k-th site whose element 
is given by 



5fc (*,0 = (-i) 5 -V* +i '- 5)(s+1)/2 i 



" 5 " 




' S ' 


i 




i' 



--: hw\S;i' - i)y 



(0 <i,i' < 5). 



[5-» + t']![5 + *-»']! \S;i' -i) k 

(3.11) 



The symbol * for two (5 + 1) x (5 + 1) vector- valued matrices 

/ boo) ••• \yos) 



is defined by 



x-ky 



( |xoo) ••• \xos) 
\ \xso) ■■■ \xss) 

( Ef=o l x o«) ® |y«o) 

V Ef=o l 3 ^} ® |yuo) 



\ 1 2/50 ) • • • \yss) 

Ef=0 l x 0n> ® \y u s) 



which is associative. 

For example, for 5 = 2, 

(^oo|2;0)fc h i\2;l) k h Q 2\2;2) k 

h 10 \2;-l) k /iii|2;0} fc /ii 2 |2;l) fc 

/i 20 |2;-2) fc /iai|2;-l) fc h 22 \2;0) k 

and the product in the form (|3.10|) is calculated as 
9i * • • • * 9L 

G/l *---*<?l) (0,0) ( 5l *---* 5L )(0,l) (^i*...*5x)(0,2) 

(gi*---*gL) (1,0) (gi*---*g L ) (1,1) (#i *•••*#£) (1,2) 

(<7i*- ■■*<?£) (2,0) (^★•••★p L )(2,l) (^★■■■★ 5l )(2,2) 



with 



(5l*"- *Sl) (*,0 

^ ] hii 2 hi 2 i 3 ■ ■ ■ h 
i fc =0,l,2 

x |2; « 2 - *)i ® I 2 ; *3 - ^2)2 ® ■ ■ ■ ® |2; ii, - i L -i)L-i ® |2; £' - 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



Then the matrix product ground state (|3.10p is 

(Si*- (0,0) + (0i*---*0 L )(l,l) + (si*..-*0£)(2,2) 

— ^ ^ hi 1 i 2 hi : 



ife=0,l,2 



(3.17) 



x |2;i 2 - n)i ® |2;i 3 - i 2 ) 2 



We define g£ by replacing each ket vector in the matrix by its corresponding bra 
vector: 



|2; - <8> |2; n - 



For example, for 5 = 2, 



^oofe(2;0| /ioifc(2;l| ^02 fc (2; 2| 
/»io fc (2;-l| h nk (2;0\ h 12 k {2;l\ 
^20fc(2;— 2| h 2 i k{2;— 1 ^22 fc (2; 0| 



(3.18) 



(3.19) 



Now we introduce "G matrix" , which will play an important role in our study. Let us set 
an (S + l) 2 dimensional vector space W and its dual orthogonal space W* as 



W= C\a,b)), W*= C((a,b\. 



(3.20) 



Here, {\a,b}} \ a,b = 0, . . . , S} ({{{a,b\ | a, b = 0, ...,S}) is an orthonormal (dual or- 
thonormal) basis. We define an (S + l) 2 x (S + l) 2 matrix G acting on the space W 
as 

G { aMc,d) = b \ G \ c ' d )) = 9 j (a, c)g(b, d), (3.21) 

or equivalently as 

G = g ] g. (3.22) 
We also introduce Ga for an operator A acting on the one-site vector space Vs as 

(GA) (aMc4) = ({a,b\G A \c,d)) = gHa,c)Ag(b,d). (3.23) 
Each element of the matrix G can be expressed explicitly as 



G(a,b;c,d) ~ $c-a,d-bT a bcd, 



(3.24) 



where 

T abcd = h ac h M = ^ ir+bq (a +b+ c +d -2S)(S + l)/2 



X 



S 
a 



S 
b 



s 

c 



s 
d 



[S - a + c)\[S + a - c]\[S - b + d]\[S + b - d]\ 



(3.25) 
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Each element of Ga for A = S z , S + and S , which act on IS; m) as 

S z \S;m) =m\S;m), 



S+\S;m) =y/(S -m)(S + m + l)\S;m + l), 
S~ \S; m) =\/ (S + m)(S — m + 1) \S\ m — 1), 

can be also expressed as 

( G S^)( atb - C: d) =fic-a,d-b(d ~ b)T abcd , 



(Gs + ) {a , b ;c,d) =$c- a ,d-b+iV(S -d + b)(S + d-b+l)T, 



(3.26) 
(3.27) 
(3.28) 

(3.29) 
(3.30) 
(3.31) 



abed ! 

( G S- )(a,b;c,d) = S c-a,d-b-l \/( S + d ~ b )i S ~ d + b + l)T abcd . 

The squared norm of the ground state is calculated as 

=Tr [»J ★■■■★»{,] Tr[ 9l *---*g L ] 
= Tr (g\*---*g\j <8> (5i*-"*5l) 
= Tr ® 01) * • • • * ® 5i) 
= TrG L 

Note that the elements of gt <8> gfc = G are no longer vectors, and thus we can replace the 
symbol * by the usual product in the third line of (|3.32p . The one point function (A) of 
an operator A can be written in terms of G and Ga as 



(3.32) 



(A) 



Tr[A 1 g 1 *g 2 *---*gL] TrG j4 G i " 1 



(*l*> 



Tr G L 



(3.33) 



where Af-gu is defined by (Akgk) (*>*') = -Afc (fffc(M'))- In the same way, the two point 
function of ^4 and £? can also be written in terms of G, Ga and Gb as 



(Tr G^) ^Tr G A G r - 2 G B G 



iL-r 



(3.34) 



Investigating the eigenvalues and eigenvectors of the matrix G will be crucial for the 
analysis of correlation functions. In the next section, we study the G matrix in detail. 

4 Spectral structure of the G matrix 

In Ref. [23], we conjectured that the spectrum of G is given by 



x e = (-iY my 



2S + 1 

S-l 



(£ = 0,1,..., S), 



where the degree of the degeneracy of each Xi is 2£+ 1. One can easily find that 

|A | > |Ai| > • • • > |A S |. 



(4.1) 



(4.2) 
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In this section, we prove the conjecture by giving an exact form for the eigenvector corre- 
sponding to each eigenvalue. 

First one observes that the G matrix has the following block diagonal structure: 

G= G {j \ G^ G EndWj, (4.3) 

-S<j<S 

C\i,i + j)) j>0, 



w= Wj, w j = < 

-s<j<s 



0<i<S-j 

C\i-j,i)) j<0. 

{0<i<S+j 



(4.4) 



The size of each block G® is (5 - \j\ + 1) x (S - \j\ + 1). Each element of is 
+ j\G®\c,c + j)) =(-l)V a+c+i_5)(S+1) [S - a + c]\ [S + a - c]\ 



(a, a 





' s ' 




S 




' S ' 




S 




a 




a + j_ 




c 




c + j _ 



(4.5) 



We construct intertwiners among the 2S + 1 blocks G^ {j = — S, . . . , S). This helps us 
to construct eigenvectors of each block from another block with a smaller size. (The same 
idea was used in Ref. [28] to study the spectrum of a multi-species exclusion process). Let 
us define a family of linear operators {Ij}-s<j<-i,l<j<S as 

Ij € Hom(^,^_i), 



((a,a + j - l\Ij\c,c + j) 



a [a + j] [S-a-j + 1] 

V ms-j+i] 



[a][S-a + l] 

q V bi - i + 1] 



a, 



a — 1, 



otherwise 



(4.6) 



(4.7) 



for 1 < j < S, and 



Ij G Bom(W j ,W j+1 ), 



f„-a [a-j][S-a + ; + l] 



((a - j - l,o|J,-|c-j,c}) 



[-i] [s + j + i] 



A- a+j 



' [a][S-a + 1] 
[-J] [S + j + l] 



c = a, 



c = a — 1, 



otherwise. 



(4.8) 



(4.9) 



for —5 < j < —1. By direct calculation, one finds 

Proposition 4.1. TTie matrix Ij enjoys the intertwining relation 

IjG® =G ( - j - 1) I j for l<j<S, 
LG (j) =G {j+1) L for - S < j < -1. 



(4.10) 
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With the use of Proposition 14.11 one can show the following: 
Theorem 4.2. Each block G^' has a simple (nondegenerated) spectrum 

Spec = {\t}\j\<e<s, 
and the corresponding eigenvectors are given by 



(4.11) 



o<i<s-e 



E * ( 

0<i<S-£ 





(\ 


[i + £]\[S-i]\ 


[S)\ 


M 


[S-i-£]\\i]\ 


l[s- 


(} 


[i + £]\ [S-i]\ 



[ S }\[£]\[s-i-e]\[{\ 



i,i + £)) j>0, 



i + e, i)) j < o, 



(4.12) 



for i = \j\, and 
for \j\ + 1 < I < S. 



Ij+i\Xe)) j+1 = Ij+iIj+2 ■ ■ ■ It\*i)) t 3 ^ °' 
= I3-1I3-2 ■ ■ ■ I-t\*t))-t J < 0. 



(4.13) 



Figure [2] is helpful to understand how the eigenvectors are constructed. We prove this 
theorem below for only j > since one can show it for j < in the same way. 

Proof of Theorem First, by direct calculation given below, we find that has 
an eigenvalue Xj and its eigenvector is \Xj))j defined by ()4.12p . Each element of G^\Xj))j 
is calculated as 

({a,a + j\G^\X j )) j 

= (-l) j q (a+C+j ~ S){S+1) [S -a + c]\[S + a-c}\ 

0<c<S-j 



S 
a 



S 
a+j 



S 

c 



-iy q 

[s]\ ""' + 

o<c<s-j 



s 




l[S-j]\[c + j]\[S-c}\ 




qU+1)C \J 


[S]\[j]\[S-c-j]\[c}\ 


' s ' 




S 


[s-j}\ 


a 




a + j 





(4.14) 



l)(5+i+2) [S-a + c]\[S + a-c]\ 
[S-c-j]\[c]\ 



Using the formula 

E 



Q<fe<n 



a + n — k 
n — k 



P + k 
k 



Ma+0+2) 



a + (3 + n + l 
n 



(4.15) 
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|As-i) 


>1 


A5-2) 


>x ••• 


|A 2 ))! IAi))! 


Ui 




I 




I 






w 


|A5» 
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|As-i) 
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1 
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|As»_ 2 


IA5-1) 
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)-2 ■■■ 


|A 2 ))_ 2 


J^-S+2 


|As))_ 5+2 


|As-i) 
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IA5-2) 


)-S+2 
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|As))_ 5+ i 


IA5-1) 


)-S+l 








t/- s 
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|A 5 ))_ 5 













Figure 2: Structure of the eigenvectors of G (|3.2ip . 



we obtain 



<<a,a + j|G^)| Aj )) 



X [5']! (? -S'j+a(j-5)+S 2 +2S+2 



(-1)^([5]!) 
Aj((a,a + j\X j )) j 



2S + 1 

s-j 




7 (i+i)» 



[S]!0![S 



Note that the first element of \Xj))j is 1 by the definition (|4.12p : ((0,j\X 



■3 ll 3 



(4.16) 



Next, we show by induction that G^> has eigenvalues Xi (j < £ < S) and their corre- 
sponding eigenvectors are given by |A^)) • defined by (|4.13p . Suppose the theorem is true for 
\Xe))j +1 ,£ = j + 1, . . . , S (j > 0), that is to say that the block has the eigenvalues A^ 

and their corresponding eigenvectors |A^)) ■ +1 [G^ +l '\X(,)).- +1 = A^|A^)) - +1 with |A^))- +1 7^ 
0) for £ = j + 1, . . . , S. Additionally, suppose that the first element of each |A^)) - +1 is 1. 



11 



Using the intertwining relation (|4.10p . one finds G^ Ij + i\\^f) - +l = \^Ij + \\\^f) We also 



/j+v 



find that the first element of Jj+i|A^))j +1 is 1, and thus Ij + i\\ifj - +1 is nonzero. Further- 
more, thanks to l\ / £2 => 7^ X^ 2 , the vectors Ij + i|A£)) J+1 + 1 < £ < S) are distinct 
(in other words, Ij+i is injective). We have already constructed the remaining eigenvector 
of G^> explicitly, which is with its eigenvalue Xj distinct from \z (3 + 1 < £ < S). □ 
The conjecture for the eigenvalues of the G matrix that we exhibited in the beginning 
of this section follows as a simple corollary of Theorem 14.21 Moreover, we constructed 
their eigenvectors which are important for computing spin-spin correlation functions. 



Proposition 4.3. The squared norm of \\ifj „■ is 



s(\j\+i)~e(e+i) 



n [ s + e + i]\[e-\j\]\[s-e\\[\j\] ] - 

[S]l[£ + \j\]\[S-\j\}\[2£ + l] 



(4.17) 



where we denote the transpose of\Ai))- by j((Xi\. 



We prove this proposition only for j > 0. 

Proof of Proposition \4-3[ One can easily show that the product of intertwiners (which 
is also an intertwiner) has the following form by induction: 



((a, a + j\Ij +1 I j+2 ■ ■ ■ h+ih\c, c + £)) 
--{-l) a - c q cj - a 



J 

■ c 



\j]\[S-t\W ] -[S-c]\[c + e\\[S-(a + 3)V. 
mS -j]\[c]\[S-a]\[a + j]\[S-{c + £)]]' 



(4.18) 



Then, j((Xe\Xt)) ,■ = e((\e\ (Ij+i •••It) Ij+i ■ ■ ■ Ii\h))t is calculated as 



[£]l) 2 [S 



3 1 ■ 0<a<S-j 
0<i,i'<S 



^ (_i) i + i 'g(*+ i ')(^+i+i)-2^ 



[S - + £}l[S - i']\[i' + £}\[a]l[S - (a + j)]l 
[i\l[S ~ (i + eWY-[S - (i' + £)}\[S -a]\[a + j\V 

The triple sum has the closed form 



[£]]) 2 [£-j]l[S + £+l]\ 
[S-£]\[j+£]\[2£ + l] ' 



which finishes the proof. 



(4.19) 



(4.20) 
□ 



5 Spin-spin correlation functions 

In the last section, we investigated the eigenvalues and eigenvectors of the G matrix. By 
utilizing Theorem 14.21 and noting (|4.2p . the one point function (^4) can be represented as 



(A) = Xo MX ^ h (-"..D 
o\\Ao|Ao); 
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in the thermodynamic limit L — > oo. As an application, we can calculate the probability 
of finding S z = m value as 



Prob(5^ =m) =(\S;m)(S;m\) 

_ [S + m]\[S - m]\ 
[25 + 1]! 



i=0 



m-5) 


5 




' S ' 




i — m 




i 



(5.2) 



The two point function (|3.34p can be also represented as 



£=0 



o((A |Ao}} oi ((A,|A,}) 



(5.3) 



in the thermodynamic limit. Inserting (|3T29L (|330|) . (EOTjl . flU]), (j4TT2l) . <^W\ and (ETTTj) 
into (|5.3p . one finds the large-distance (r — > oo) behaviors of the spin-spin correlation 
functions (SfS*) and (S^S~) are 



[3] [S + 2] 



g 25 - 2 [5]([25 + l]!) 2 

[2] [3] [5 + 2] 
g 35 - 2 ([25 + l]![5]) 2 



(o((Ai|G 5 ,|Ao)} ) 2 
U<<Ai|G s -|Ao>>o) 2 



<sysjf) 

where 

o((Ai|G s .|Ao» 

x {g 5+1 + g- 5 - 1 - (g + g _1 )9 2i '- ,s }[S + i - + i' - *]! 



[5 + 2] 
[51 



[5 + 2] 



(5.4) 
(5.5) 



? -s 2 -s-i 



i,i'=0 
-l\„2i'-Si 



/\„(S+2)(i+i') 



(5.6) 





" 5 " 




" 5 " 


*]! 






i' 



-i((A 1 |G 5 -|Ao»o = o«Ao|G s+ |Ai))_ 1 

5 5-1 

= _ ? -S 2 -S/2+l/2 ^ ^ (? (S+2)i+(S+3)i' ^ 
i=0 i'=Q 

x ^/(5 + i - + i - i'](5 - i + V + 1)[5 - i + z' + 1] 

" 5 " 



5 




" 5 " 


. *" + 1 . 




i' _ 



(5.7) 



x 



y/\i f + l][S-i'}[S]-i[S + i' - i}\[S + i - »' - 1]! 



Note that the terms with (j,£) = (0,1) and (—1,1) in (|5.3p dominate the large-distance 
behaviors of {SfS*} and (5^5~), respectively, since 

o«Ao|G<HAo)) = i«Ai|G 5 »|Ao» = -i((M\G s * |A )) = 0, (5.8) 
o((A |G s -|Ao)) = i«Ai|G 5 -|Ao» = o«Ai|G 5 -|A » = 0. (5.9) 

Both (5f S*} and (SfS~) exhibit exponential decay with correlation length 



C = In 



[5 + 2] 



-i 



(5.10) 



generalizing the results for q = 1 [20j or 5 = 1 [7] case. 
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6 Conclusion 



In this paper, we investigated one and two point functions of the q-VBS ground state 
of an integer spin model (the (/-deformed higher-spin AKLT model). The formulation 
of correlation functions by use of the matrix product representation of the ground state 
shows that the structure of a matrix, which we call G matrix, plays an important role. We 
obtained the eigenvalues and eigenvectors of the G matrix with the help of constructing 
intertwiners connecting different blocks of G. Then we calculated the spin-spin correlation 
functions by use of the eigenvalues and eigenvectors of the G matrix, and determined the 
correlation amplitudes and correlation lengths of the longitudinal and transverse spin-spin 
correlation functions. 

It is interesting to investigate other types of correlation functions. For example, the 
entanglement entropy, which is defined in terms of the reduced density matrix, is a typical 
quantification of the entanglement of quantum systems. It is intriguing to calculate the 
entanglement entropy for the g-deformed model and observe the change from the isotropic 
point I2DE21I3D1 (see also Refs. [SUES] for other VBS states). 



Acknowledgements 

The authors thank Atsuo Kuniba and Kazumitsu Sakai for useful discussion. CA also 
thanks Kirone Mallick and Andreas Schadschneider for the very kind hospitality during his 
stay in Europe. This work is supported by Grant-in- Aid for Young Scientists (B) 22740106 
and Global COE program "Education and Research Hub for Math-for-Industry." 



A Proof of Proposition 12.1 



We prove Proposition|24Jby induction. Suppose that Proposition|2_ L l]holds for ^AA^ I / /. 



We calculate the four terms of the action of 

AXZa = 



q-q~ l x a q l/V? 



-1 rr q- 1 y/q 1/y/q 



q-q 1 x c 



(A.1) 



on (^AX a/3 S j vj, separately. 



^D*"®D x iD^, J (AXZ.Y vj I ' q^ s \n]\{x a xp) J ^ n+1 ^ 



2//5+J-n 


" J ' 




J 




. V _ 




n — fj, 



2S-J 



2(v-S-l) 



(A.2) 



XpVa 
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(g^ VlfJpiQ (^)\j/ q ^nn)K*c*0) J -™ (A-3) 

25- J 



At=0 



2^5-2/^-45+ J+n 


' J ' 




J 




. V _ 




n — u 



. ^2 q-2/J,S-2n-2S+J+n 





J " 




J 


A* _ 




n — 



(A.4) 



2S-J 



D*°? ,-D 



^®^x)(A^)% J /^ 5 [n]!(x a x / ,) J -("+ 1 ) 



(A.5) 



2/t5-25-J+n 


" J " 




J 








n — /i 



E< 

IPC3|I - (TOD gives 



E< 

/i=0 



n 

E* 

At=0 



2S-J 



{x a ypY +l (xpy a ) n M ] [ (x a y p - q 



_2(i/-S-l) 



-2nS-2pL-2S+J+n 


" J " 




J 




. a* . 




n — fi 




25- 


-J 


x a yp - q 2S xpy a ) ] | 


[x a yp - 




V- 


=i 




-2/j.S- 2^-25+ J+n 


' J ' 




J 




. A 4 . 




n — u 



,2(>-5) 



X0c 



{x a ypY (xpy a ) 



n—fi 



2S-J 



[x a yp-q 2{S J) X0 a \ \ [ [Xayp-q 2 ^ S 1) x p y c 



,14=0 



2/i5-2/i-25+J+n 


" J " 




J 




. _ 




n — a 



{x a ypY +l {X0a) 



n— n 



2S~J 



H=0 

25- J 



2(^-5-1) 



2fj,S~2fi-J+n 


' J ' 




J 




. A 1 _ 




n — u 



i 



(x a ypY (xpy a ) 



n+1— /i 



(A.6) 



i/=i 
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Dividing the first term of ()A.6|) — (|A.5j) by q — q 1 , we obtain 

j2 q -^s-,-2s +n[J _ M] r j j r j i (wr „ 

LX U LI 

fj,=0 J 

2S-J 

x IT ( x <xy&-i {v ~ s ~ x) ^0a 



n+l 

: ^ 9 -2„S- M +n+l[j _ 
i/=l 



J 




J 


. M-l _ 




n. + 1 — fj, 



(A.7) 



where we replaced // — > /i — 1. Dividing (]A.2[) — the second term of (|A.6|) by g — g 1 , we 
obtain 





" J " 




J 


n] 










. V _ 




n — \x 



{xaypY ( x py 



H i \n+l—(i 



2S-J 



(A.8) 



i/=l 



Since 



J 



0, we can extend ^ in ()A.7|) and (|A.8|) to < u < n + l. Finally we have 



n+l 



453 + (Ja^D = q~ 2flS [™ + 1] 

fj,=0 

2S-J 



' J ' 




J 


. /" _ 




n + 1 — /i 



(A.9) 



IT \ x <±yp - i 2{v s 1)x 0a 



and since 



(AX^y + \j j q^ S [n]\(x Q x^ J -^ 



COD + ggj) - (pop - (|A3D 

da + m> - mi 



(A.10) 



(ED, 



Proposition 12.11 is true for ( AX 



n+l 
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